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ABSTRACT 


The  free  eleetron  laser  (PEL)  is  among  the  latest  technologies  of  interest  to  the 
U.S.  military,  in  particular,  the  Navy.  In  naval  applications,  PEL  laser  would  serve  as  a 
self-defense  weapon  system,  protecting  the  ship  from  an  array  of  threats  including  anti¬ 
surface  cruise  missiles  and  small  boats.  This  system’s  potential  range  and  deep  magazine 
makes  it  ideal  as  point  defense  against  incoming  missiles.  Its  inexpensive  cost  of  only  a 
few  dollars  per  engagement  and  multi-mission  capability  makes  this  future  weapon 
system  superior  to  the  short-range  missile-defense  systems  employed  today.  The  most 
powerful  PEL  is  currently  located  in  Jefferson  Lab,  operating  at  10  kW,  two  orders  of 
magnitude  short  of  the  1  MW  power  level  required  for  weapons  application.  This  thesis 
will  describe  the  components  and  theory  of  operation  of  the  PEL,  as  well  as  analyze  two 
competing  designs  for  the  next  step  in  the  evolution  of  the  future  weapon  system,  the  100 
kW  PEL,  proposed  by  Brookhaven  and  Los  Alamos  National  Labs.  Due  to  advances  in 
NFS  simulation  techniques  for  the  amplifier  configuration,  a  more  in  depth  analysis 
including  the  effects  of  electron  beam  tilt  and  shift  is  performed  for  the  first  time  on  these 
proposed  designs. 
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I.  INTRODUCTION 

A.  HISTORY  OF  THE  FREE  ELECTRON  LASER 

The  free  electron  laser  (FEL)  is  among  the  latest  technologies  of  interest  to  the 
U.S.  military,  in  particular,  the  Navy.  In  naval  applications,  FEL  laser  would  serve  as  a 
self-defense  weapon  system,  protecting  the  ship  from  an  array  of  threats  including  anti¬ 
surface  cruise  missiles  and  small  boats.  This  system’s  potential  range  of  six  to  ten 
kilometers  and  deep  magazine  makes  it  ideal  as  point  defense  against  incoming  missiles. 
Its  inexpensive  cost  of  only  a  few  dollars  per  engagement  along  with  a  multi-mission 
capability  versus  using  a  multi-million  dollar  missile  per  engagement  makes  this  future 
weapon  system  superior  to  the  short-range  missile-defense  systems  employed  today. 

The  FEL  is  another  product  in  a  long  line  of  methods  of  generating  coherent 
electromagnetic  energy,  all  of  which  prove  quite  useful  today.  The  first  of  these 
techniques,  developed  in  the  1880’s,  was  the  antenna.  The  antenna  emits  energy  in  the 
form  of  radio  waves  directionally  and  coherently,  generated  by  oscillating  currents  in 
metal  and  wire.  The  primary  purpose  of  the  antenna  was  the  transmission  and  receipt  of 
radio  signals.  Through  the  years,  antennae  have  retained  this  function  in  communications 
and  have  been  utilized  in  the  entertainment  industry  in  the  service  of  radio  and  television. 
They  now  serve  as  the  backbone  of  all  wireless  telephone  and  broadcast  television 
networks  worldwide. 

The  concept  of  coherent  emission  was  later  applied  to  the  microwave  region  of 
the  electromagnetic  spectrum  in  the  creation  and  development  of  microwave  tubes.  The 
microwave  tube  was  first  used  in  the  1930’s,  primarily  for  communication  and  radar 
systems.  The  microwave  tube  marked  the  beginning  of  coherent  energy  generation  in  a 
closed  container.  The  radiation  in  the  closed  cavity  interacted  with  the  relatively  slow 
moving,  non-relativistic  electrons.  The  resulting  radiation  had  long  wavelengths  on  the 
order  of  millimeters.  The  tubes  were  very  efficient,  utilizing  better  than  half  of  the 
energy  supplied  as  input.  The  most  common  use  for  these  tubes  now  is  in  an  appliance 
found  in  the  kitchens  of  almost  every  home  nationwide,  the  microwave  oven.  Building 
and  improving  upon  the  concept  of  microwave  tubes,  the  Microwave  Amplification  of 
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Stimulated  Emission  of  Radiation,  or  MASER  was  built  and  developed  by  Charles  H. 
Townes  in  1954  [1], 

The  stimulated  amplification  of  radiation  was  later  used  to  produce  shorter 
wavelength  radiation  in  the  visible  region  of  the  spectrum.  This  gave  rise  to  the  Light 
Amplification  of  Stimulated  Emission  of  Radiation  or  what  we  commonly  call  the  laser. 
Theodore  H.  Maiman  built  the  first  working  laser  in  1960  [1].  Throughout  the  1960’s, 
several  types  of  lasers  were  produced  using  electrons  bound  in  various  atoms  and 
molecules  such  as  ruby,  uranium,  helium-neon,  semiconductor  materials,  and  carbon 
dioxide.  The  conventional  laser  now  comes  in  all  shapes  and  sizes  and  is  commonplace 
in  society,  including  laser  pointers,  printers,  compact-disc  players,  DVD  players,  product 
scanners  and  even  as  instruments  in  medical  procedures. 

With  all  the  promise  of  the  conventional  laser,  it  also  had  shortcomings.  These 
lasers  had  short  wavelengths,  but  they  were  not  tunable.  Bound  electrons  transitioning 
from  excited  to  lower  energy  states  determined  the  wavelengths  of  these  lasers.  They 
were  restricted  by  the  energy  levels  of  the  atom  or  molecule  that  the  laser  is  based  upon. 
Adjustment  of  wavelength  was  not  possible.  Furthermore,  these  lasers  were  not  as 
efficient  as  their  predecessors  in  the  microwave  regime. 

The  issue  of  tunability  and  efficiency  of  lasers  was  addressed  in  1971  by  John 
Madey  [2].  His  invention,  the  FEE,  combined  the  benefits  of  both  the  earlier  lasers  and 
the  MASER,  providing  short  wavelengths  and  relatively  high  efficiency.  With  the 
evolution  of  technology  through  the  years,  the  free  electron  laser  has  been  improved  upon 
throughout  the  world  and  has  most  recently  achieved  a  power  level  as  high  as  ten 
kilowatts.  This  laser  is  currently  in  operation  at  Jefferson  Lab  in  Newport  News, 
Virginia. 

B.  ADVANTAGES  OF  THE  FREE  ELECTRON  LASER 

One  of  the  major  advantages  of  the  FEL  over  other  lasers  is  its  flexibility  of 
design.  In  conventional  lasers,  the  operating  wavelength  is  predetermined  by  the 
chemical  composition  of  the  lasing  medium.  In  contrast,  the  basic  components  of  the 
FEL  can  be  used  in  the  microwave,  infrared,  visible,  ultraviolet,  and  X-ray  regions  of  the 
electromagnetic  spectrum.  Furthermore,  with  a  given  configuration,  the  free  electron 
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laser  has  an  ability  to  be  finely  tuned  in  operating  wavelength.  To  accomplish  this,  one 
would  only  have  to  change  the  electron  beam  energy  through  adjustment  of  the  intensity 
of  the  electric  fields  in  the  radio  frequency  (RF)  accelerator  or  injector.  Free  electron 
lasers  are  capable  of  running  twenty-four  hours  a  day  for  weeks  at  a  time.  Another 
advantage  of  the  free  electron  laser  over  conventional  lasers  is  cost  of  operation.  For 
example,  in  high-power  conventional  chemical  lasers,  the  chemical  being  used  as  the 
lasing  medium  is  spent  and  exhausted  during  operation.  Since  the  chemical,  usually 
toxic,  is  being  expelled,  it  must  be  replenished  after  repeated  use  for  the  laser  to  operate. 
The  FEL  does  not  suffer  from  this  shortcoming.  As  long  as  sufficient  power  is  supplied 
to  the  system,  and  none  of  the  components  are  damaged,  the  free  electron  laser  will 
operate.  There  is  no  exhaust  gas  involved  in  the  process  of  lasing,  thus  no  need  for 
replenishment  and  no  danger  to  people  or  the  environment. 

The  one  drawback  of  the  FEL  is  the  fact  that  the  system  tends  to  be  quite  large 
and  expensive,  but  only  compared  to  the  conventional  laser  at  low  power.  Any  MW- 
level  conventional  laser  will  have  a  comparable  large  cost.  The  size  of  the  FEL  does  not 
scale  linearly  with  output  power  (i.e.  a  lOkW  FEL  is  not  one-tenth  the  size  of  a  lOOkW 
FEL).  Furthermore,  advances  in  FEL  designs  incorporating  technologies  such  as 
superconducting  RF  accelerators  and  application  of  concepts  including  short  Rayleigh 
lengths,  will  reduce  the  size  of  the  system. 
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II.  FREE  ELECTRON  LASER  SYSTEM  COMPONENTS 


A.  INTRODUCTION 

There  are  several  components  that  make  up  the  total  free  electron  laser  system. 
Many  supporting  elements  are  used  for  cooling,  input  power,  and  other  functions,  but  for 
brevity  the  focus  will  be  limited  to  the  major  components  of  the  system.  The  basic 
system  is  comprised  of  an  injector  to  provide  the  electron  beam,  an  accelerator  to  raise 
the  electron  beam  to  a  useful  energy  level,  an  undulator  to  extract  light  energy  from  the 
electron  beam,  a  resonator  to  store  light  created,  a  beam  dump  to  dispose  of  the  used 
electrons,  and  bending  magnets  to  route  the  electron  beam  through  the  system  to  each  of 
components.  A  schematic  of  the  free  electron  laser  system  at  Jefferson  Laboratory  is 
shown  in  Figure  1  below. 


UV  wiggler 


Figure  1 .  Diagram  of  Jefferson  Lab  free  electron  laser  system  ([From  Ref  [3].) 

B.  INJECTOR 

The  electron  beam  formation  begins  at  the  injector  where  free  electrons  are 
emitted  from  a  cathode.  Four  types  of  injectors  are:  thermionic,  direct  current  (DC), 
radio  frequency  (RF)  and  superconducting  RF.  The  thermionic  gun  utilizes  heat  energy 
to  excite  electrons  in  a  cathode  (a  metal  plate),  providing  the  source  of  electrons  used  in 
the  system.  This  cathode,  when  heated,  releases  electrons.  The  RF  and  DC  injectors 
accomplish  the  same  using  the  photoelectric  effect;  a  laser  illuminates  the  cathode,  and 
the  cathode  absorbs  the  light.  This  absorbed  energy  excites  electrons  in  the  cathode, 
causing  it  to  release  the  electrons.  The  newly-freed  low  energy  electrons  are  given 
additional  energy  and  accelerated  to  relativistic  velocities  by  the  injector.  The  beam  is 
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typically  accelerated  to  only  modest  energy  levels  (2-10  MeV)  in  the  injector.  A 
schematic  of  the  injector  is  shown  below  in  Figure  2. 
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Figure  2. 


7-10  MeV 

Schematic  for  an  RF  injector  gun  ([From  Ref  [4].) 


C.  ACCELERATOR 

After  leaving  the  injector,  the  free  electrons  are  relativistic  but  still  remain  at  too 
low  an  energy  to  be  useful  for  weapons  application.  The  accelerator  (shown  in  Figure  3 
below)  gives  the  electron  pulses  the  additional  energy  needed  for  the  desired  power 
output  of  the  laser.  The  accelerator  is  a  series  of  cavities  into  which  RF  energy  is  fed  and 
alternating  electric  fields  are  set  up.  The  size  of  cavities  roughly  determines  the 
frequency  of  operation  of  the  system.  The  number  of  cavities  in  the  accelerator  roughly 
determines  the  amount  of  RF  energy  in  the  accelerator.  During  operation,  electrons  go 
through  each  cavity  in  the  accelerator  and  are  exposed  to  the  alternating  electric  fields 
that  are  present.  When  the  electrons  are  coming  from  the  injector  in  phase  with  the 
electric  field,  they  absorb  energy  from  it.  This  energy  absorption  accelerates  and  bunches 
the  electrons,  and  slightly  reduces  the  electric  field  present  in  the  cavity,  so  the  RF  energy 
must  be  replenished  to  continue  operation.  The  power  requirement  to  replenish  the 
electric  fields  is  reduced  by  the  recirculation  of  electrons  through  the  system.  Electrons 
that  have  gone  through  the  undulator  can  be  sent  through  the  accelerator  for  a  second 
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time,  now  out  of  phase  with  the  electric  field  in  each  cavity.  The  electrons  then  lose 
energy  and  replenish  the  RF  electric  fields. 


Figure  3.  An  example  of  a  RF  accelerator  ([From  Ref  [5].) 


D.  UNDULATOR 

Once  emerging  from  the  accelerator,  the  electron  beam  is  at  a  high  energy 
(approximately  100  MeV).  The  electron  beam  energy  is  the  source  of  energy  for  the  laser 
light.  Extraction  of  energy  from  the  electron  beam  and  amplification  of  light  are 
accomplished  in  the  undulator.  The  undulator  (shown  in  Figure  4)  is  constructed  from 
two  rows  of  adjacent  permanent  magnets  with  alternating  poles.  The  magnets  are 
oriented  so  that  there  is  an  alternating  static  magnetic  field  along  the  undulator  axis. 
Upon  entry  into  the  undulator,  the  electrons  experience  Lorentz  forces  due  to  the 
magnetic  fields.  The  Lorentz  forces  cause  the  electrons  to  oscillate  with  small  amplitude. 
These  small  oscillations  cause  the  electrons  to  emit  radiation  in  the  form  of  light.  The 
light  emitted  by  the  electrons  interacts  with  the  light  stored  in  the  undulator  through 
stimulated  emission,  which  amplifies  the  light  already  present  in  the  undulator.  The 
undulator  typically  extracts  only  a  few  percent  of  the  energy  from  the  electron  beam  for 
output  laser  light.  The  rest  of  the  electron  beam’s  energy  is  recirculated  and  routed  back 
to  the  accelerator  for  increased  system  efficiency. 


7 


Figure  4.  A  schematic  of  an  undulator  ([From  Ref.  [6].) 

E.  OSCILLATOR  RESONATOR 

The  oscillator  configuration  is  equipped  with  a  resonator  to  store  and  amplify 
light  through  many  passes.  The  resonator  consists  of  an  evacuated  cavity  terminated  by 
two  specially  constructed  mirrors.  One  mirror  is  almost  perfectly  reflective,  while  the 
second,  the  out-coupling  mirror,  is  partially  transmissive,  allowing  a  fraction  of  the  light 
to  pass  through  it  and  escape  the  resonator  as  a  usable  laser  light  source.  The  undulator  is 
typically  placed  halfway  between  the  mirrors.  Light  bounces  between  the  two  mirrors  at 
the  operating  wavelength  and  is  amplified  between  each  bounce  by  the  undulator  through 
stimulated  emission  due  to  its  interaction  with  the  electron  beam. 

F.  AMPLIFIER  SEED  LASER 

The  amplifier  configuration  of  the  FEL  uses  a  long  undulator  to  amplify  a 
coherent  light  source  in  a  single  pass.  Therefore,  the  amplifier  design  has  no  resonator 
mirrors  and  no  ability  to  store  light.  The  amplifier  employs  a  seed  laser  as  a  source  of 
light  to  be  amplified  by  the  undulator.  A  low-power,  solid-state  laser  usually  serves  as 
the  seed  laser  for  an  amplifier  FEL  system. 

G.  BENDING  MAGNETS 

In  the  recirculating  electron  beam  design  of  the  FEL  system,  there  are  bends  in 
the  system  where  the  electron  beam  must  be  routed  around  components.  A  high  power 
electron  beam  incident  on  either  of  the  resonator  mirrors,  the  injector,  or  the  system 
piping  itself  would  prove  disastrous  to  the  operation  of  the  system.  Bending  magnets 
provide  the  means  of  routing  the  electron  beam  through  the  system  piping  around  these 
components.  The  powerful  bending  magnets  are  configured  so  that  the  electron  beam. 
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when  passing  through  them,  is  deflected  by  the  Lorentz  forces  resulting  from  the 
magnetic  fields.  Bending  magnets  provide  an  extremely  tight  control  of  the  electron 
beam  direction.  However,  the  full  recirculation  through  the  bending  magnets  requires 
that  the  electron  beam  has  a  small  energy  spread  to  ensure  efficient  operation. 

E.  BEAM  DUMP 

Once  the  electron  beam  has  passed  through  the  undulator  and  has  lost  energy,  it  is 
routed  through  bending  magnets  back  to  the  accelerator.  Even  though  the  electrons  have 
lost  a  few  percent  of  their  energy  in  the  undulator,  they  still  have  much  more  energy  than 
those  emerging  from  the  injector.  The  higher  energy  electrons  by  design  arrive  at  the 
accelerator  out  of  phase  with  the  electric  fields.  Going  through  the  accelerator  out  of 
phase  has  the  opposite  effect  on  the  electrons  causing  them  to  lose  energy  to  the 
accelerator.  The  accelerator  decelerates  the  spent  electrons  in  preparation  for  disposal 
bringing  them  to  a  level  roughly  that  of  electrons  exiting  the  injector.  Though  at  a  much 
lower  energy,  the  electrons  still  have  energy  that  must  be  absorbed,  which  is  the  purpose 
of  the  beam  dump.  It  is  basically  a  block  of  metal,  which  just  absorbs  electrons  incident 
on  it.  Along  with  absorbing  energy,  it  absorbs  heat  and  potentially  harmful  radiation.  A 
beam  dump  used  by  the  PEL  system  at  Jefferson  Lab  is  shown  below  in  Figure  5.  Care 
must  be  taken  to  cool  and  shield  this  component  in  the  system. 


Figure  5.  Jefferson  Lab  Beam  Dump  ([From  Ref  [4].) 
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III.  FREE  ELECTRON  LASER  BASIC  THEORY 


A.  INTRODUCTION 

Of  all  components  in  the  free  electron  laser,  the  most  physics  concepts  emerge 
from  the  study  of  electrons  and  light  propagating  together  through  the  undulator.  The 
interaction  of  the  laser  light  and  electron  beam  present  in  the  undulator  is  crucial  to  the 
proper  operation  of  an  PEL.  This  section  mathematically  describes  the  process  of  laser 
light  amplification  and  propagation  along  the  undulator. 

B.  LORENTZ  FORCE  EQUATIONS 

When  the  electrons  enter  the  undulator,  they  are  exposed  to  a  transverse,  periodic 
magnetic  field.  For  a  helical  undulator,  this  magnetic  field  is, 

B  =  5(cos(A:qz),  sin(A:oz),  0) .  (HI.  1) 

In  this  equation,  B  represents  the  magnetic  field  strength  and  is  the  wavenumber 
corresponding  to  the  undulator  period  =271 1 X^).  The  undulator  period  is  the 
distance  between  two  successive,  identically  oriented  magnets  in  the  undulator. 

Along  with  the  magnetic  field  due  to  the  undulator  magnets,  electrons  are 
influenced  by  the  magnetic  and  electric  fields  of  the  laser  light,  which  (in  cgs  units)  are 
represented  by 

E  =  £'(cos  If/,  -  sin  \f/,  0) ,  (III-2) 

5  =  £'(sin^,cos^,0),  (III. 3) 

y/ =  kz  -  cot  +  (/) .  (III-4) 

In  these  equations,  E  represents  the  optical  electric  field  amplitude,  k  is  the 
wavenumber  related  to  the  optical  wavelength  X{k  =  27r ! X),  co  =  kc  =  27vcl X  is  the 
optical  frequency,  and  cj)  is  the  optical  phase. 

A  charged  particle  going  through  electric  and  magnetic  fields  encounters  a  force 
that  acts  upon  it  given  by  the  Lorentz  force  equation.  For  relativistic  particles,  the 
complete  Lorentz  FEE  force  equations  for  electrons  in  the  undulator  are  given  by 
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(IIL5) 


diyh 

dt 


me 


= 

r^=\-p\ 


(IIL6) 

(111.7) 


In  these  equations,  -e  represents  the  electron  charge,  m  the  electron  mass,  E  is  the 
electric  field,  B  is  the  magnetic  field,  P  is  dimensionless  and  proportional  to  the 
electron  velocity  =  v ! c),  where  c  is  the  speed  of  light,  and  y  is  the  relativistic 
Lorentz  factor  appearing  in  the  total  electron  energy  . 

C.  ELECTRON  MOTION  IN  THE  ABSENCE  OF  LIGHT  IN  UNDULATOR 

Electron  motion  in  the  undulator  is  first  studied  without  the  presence  of  light. 

Since  there  is  no  light  in  the  undulator,  there  is  no  electric  field  ( £  =  0 ),  and  the  force 
equations  revert  to  the  simplest  form.  With  no  light  in  the  undulator,  equations  IIL5  and 
IIL6  become 

^^A  =  -—{PxB),  (IIL8) 

dt  me 


(111.9) 


Equation  IIL9  shows  that  the  magnetic  field  does  no  work  on  the  electrons  and  the 
electron  energy  ( yme^ )  remains  constant.  Substituting  the  helical  undulator  magnetic 
field  (equation  III.l)  into  equation  III.8  results  in 

(-A  sin(V>.  A  cos(i,z),  P,  sin( V)  -  P,  cos( V»  •  (HI- 1 0) 

dt  me 

Equation  III.  10  shows  that  the  motion  of  the  electrons  in  both  the  x  and  y  directions  is 
sinusoidal  in  nature  with  only  a  phase  difference  between  them.  Integrating  the 
transverse  components  of  equation  III.  11  gives  the  transverse  velocity  ={Px^Py-‘^)) 
of  the  electrons. 
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h=-- - ^(cos(V),sin(V),0).  (III.ll) 

iTT/mc 

In  equation  III.ll,  all  integration  constants  are  zero,  indicating  perfect  injection  into  the 
undulator.  A  more  compact  form  of  this  equation  is  given  by  equation  111.12,  where 
K  =  eBAg  /  is  the  dimensionless  undulator  parameter.  The  value  of  the  undulator 

parameter  is  typically  on  the  order  of  unity  (iC  « 1 ),  so  that 

P^=-— {cos{k^z),  sin(A:oZ),  0) .  (III.  1 2) 

r 

Since  the  electrons  in  the  undulator  are  highly  relativistic,  the  velocity  of  the  electrons 
can  be  approximated  as  the  speed  of  light.  Therefore,  the  approximation 
( z{t)  «  /3ft «  ct )  is  reasonable.  Substitution  of  this  approximation  yields 

^  K 

- (cos(ryot),  sin(ryot),  0) .  (III.  13) 

r 

In  this  equation,  co^  =  kf  represents  the  electron  oscillation  frequency.  Integration  of 
equation  III.  13  describes  the  transverse  motion  of  the  electrons  in  the  undulator, 

KX 

^  (-  sin(<»ot),  cos(<»ot),  0) .  (III.  1 4) 

Itt/ 

D.  ELECTRON  MOTION  IN  THE  UNDULATOR  WITH  LASER  LIGHT 
INTERACTION 

Both  motion  and  energy  exchange  of  the  electrons  in  the  undulator  with  laser  light 
can  be  described  by  using  the  aforementioned  complete  relativistic  Lorentz  PEL  force 
equations  with  the  magnetic  and  electric  fields  due  to  the  light  in  the  undulator. 
Substitution  of  the  magnetic  and  electric  fields  into  equation  III. 5  yields 

=  — —  [£■(!  -  )(cos  -  sin  0)  +  P^B{-  sm/kf),  cosikf),  0)] ,  (III.  15) 

dt  me 

cos  yz-p  sin  y/)  +  B{p^  sm{kf)  -  p  co^ikf))] .  (III.  16) 

dt  me 
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In  these  equations,  (5^  represents  the  magnitude  of  the  longitudinal  z  component  of  the 
electron  velocity,  represents  the  transverse  x  and  y  components,  and  and  are 
the  X  and  y  components  of  the  velocity,  respectively. 

The  electrons  going  through  the  undulator  are  highly  relativistic,  so  their 
velocities  along  the  axis  of  the  undulator  z  -axis  are  very  close  to  the  speed  of  light.  This 
causes  to  be  very  close  to  but  less  than  one.  Applying  this  condition  to  equation 

III.  15  makes  the  electric  field  term  negligible  compared  to  the  magnetic  field  term. 
Integration  of  equation  III.  15  approximately  results  in  equation  III.  11  and  one  can 
proceed  as  in  the  case  of  electron  motion  in  the  absence  of  light  to  describe  the  transverse 
motion  of  electrons  in  the  undulator. 

In  studying  the  microscopic  variation  of  motion  of  the  electron  along  the 
undulator  axis,  equation  III. 6  is  utilized,  substituting  the  transverse  velocity  (equation 
III.  12)  for  and  The  result  is 

—  =  Y  = - p-E  = - (cos  k^z  cos  y/  -  sin  kf^zsmy/)  = - cos(i^  +  (p) ,  (III.  1 7) 

dt  me  yme  ymc 

where  (^  =  (k  +  k^z  -  cot  represents  the  electron  phase.  The  electron  phase  follows  the 

microscopic  position  with  respect  to  the  undulator  field  and  the  co-propagating  optical 
field  and  is  proportional  to  z(t) .  The  electron  phase  closely  corresponds  to  the  electron’s 
position  within  an  optical  wavelength.  The  first  and  second  derivatives  of  the  electron 
phase  are  taken  with  respect  to  the  dimensionless  time  T  =  ct !  L, 

dr  dt  c 

^  d^<;P 

dP  ~  dP  c"  ’ 

where  L  is  the  length  of  the  undulator.  The  time  derivative  of  the  electron  phase  gives 
the  electron  phase  velocity 


v  =  ^  =  L[{k  +  kpp^-k]. 
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An  expression  for  /?  is  found  by  squaring  the  expression  for  the  transverse  velocity 
(equation  111.13),  combining  that  result  with  the  square  of  the  z  component  of  the 
velocity  =  Ax  +  substituting  into  equation  IIL7.  The  result  is  the  following 

equation,  relating  the  electron  energy  to  its  z-component  of  velocity. 


L  (1+^^) 

i  ■ 


For  highly  relativistic  electrons  (y  » 1 ),  the  above  equation  can  be  approximated  as 


A«i- 


l  +  K^ 


(111.18) 


For  relativistic  electrons  near  resonance  A: »  ,  the  first  derivative  of  the  electron  phase 

then  becomes 


{k  +  k^) 

-k 

1  J 

^-Lk 


1  +  i^^^ 

J 


Taking  the  second  derivative  yields 


ykL 


r 


(III.  19) 


Substitution  of  equation  III.  17  for  y  in  the  above  equation  gives  the  dimensionless 
pendulum  equation. 


J  =  v  =  |a|cos(^  +  <A),  (III.20) 

where  |a|  represents  the  dimensionless  optical  field  amplitude  \a\  =  A7rNeKL\E\l y'^mc^ . 

Values  for  \a\  that  are  less  than  n  are  considered  weak  fields  and  those  that  are  greater 

than  or  equal  to  k  are  strong  fields.  This  dimensionless  pendulum  equation  describes  the 
microscopic  motion  of  electrons  within  the  undulator. 
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E.  FEE  RESONANCE  CONDITION 


Equation  III.  17  describes  not  only  the  longitudinal  motion;  it  also  describes  the 
energy  exchange  between  the  electrons  and  the  laser  light  in  the  undulator.  If  the  factor 
cos(4'  +  4>)  is  positive,  then  the  electron  energy  increases;  if  it  is  negative,  the  energy 
decreases.  For  effective  energy  exchange,  the  factor  must  not  oscillate  rapidly  in  the  time 
that  the  electron  takes  to  go  through  the  undulator,  L/PzC  ~  L/c.  Electron  phase  velocity 

{v  =  <^LI  c  =  L[{k  +  k^l3^  -  A:)]  =  )  is  the  rate  of  change  of  the  electron  phase  with  respect 

to  the  dimensionless  time. 

When  the  phase  velocity  is  equal  to  zero,  the  free  electron  laser  is  at  resonance 
and  there  is  an  optimal  energy  exchange  between  the  electrons  and  laser  optical  field.  It 
occurs  when  the  electron  velocity  along  the  axis  of  the  undulator,  fie,  is  equal  to 

kc/(k  +  k^).  For  high-energy  electrons,  J3^»l  and  so  the  corresponding 

wavelength  of  the  laser  light,  X  =  2n/k,  is  much  smaller  than  the  period  of  the  undulator, 
Xo  =  2n/ko.  In  addition,  the  oscillation  frequency  of  the  laser  light,  o)  =  Ind  X  ,  is  much 
greater  than  the  electron  oscillation  frequency,  co^  =  Inc  I in  the  undulator.  Solving 
the  resonance  condition  for  laser  wavelength  gives 

X  =  X^{\-fi:)l  (5^.  (III.21) 

A  more  useful  form  of  the  resonance  equation  is  obtained  by  substitution  of  equation 
III.  18  into  equation  III.21.  The  result  is 

X^X^{\  +  K  )  ^  (IIL22) 

2y 

F.  FEE  PHASE  SPACE 

In  a  typical  free  electron  laser,  there  are  approximately  a  million  electrons 
randomly  spread  over  each  wavelength  of  light  in  the  electron  beam.  Each  electron  in  the 
electron  beam  has  its  own  initial  phase  and  phase  velocity  (Q,v^)  when  entering  the 
undulator.  The  electron  phase  measures  the  position  of  the  electron  in  the  beam  on 
the  scale  of  the  length  of  the  laser  wavelength,  while  the  phase  velocity  is  proportional 

to  the  energy  of  the  electron.  By  examining  the  FEE  pendulum  equation  (equation 
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111.20),  one  sees  that  the  electron  will  either  gain  or  lose  energy  depending  on  its 
particular  initial  phase  and  phase  velocity.  To  graphically  aid  in  the  study  of  the  electron 
evolution  within  the  undulator,  the  electrons  are  plotted  using  a  phase  space  consisting  of 
the  electron  phase  and  phase  velocity  This  method  graphically  identifies  the 

electrons  losing  energy  (decreasing  v)  or  gaining  energy  (increasing  v),  and 
simultaneously  shows  how  the  electrons  evolve  in  phase  ^ .  Given  initial  conditions,  the 
electron  path  is  predetermined  in  phase  space  in  accordance  with  the  pendulum  equation. 
Substitution  of  the  initial  conditions  and  application  of  the  pendulum  equation  gives  the 
phase  space  path  of  an  electron 

^  ^  +  2|a|[sin(^ sin(^g  +  <j))\  .  (III.23) 

Since  some  of  the  electrons  are  gaining  energy  while  some  are  losing  energy,  electron 
bunching  in  each  wavelength  may  and  often  does  occur.  This  phenomenon  is  easily  seen 
with  the  use  of  phase  space.  Figure  6  is  an  example  of  a  phase  space  plot  of  an  idealized 
FEL  with  20  sample  electrons  exactly  at  resonance  (v,,  =0)  with  zero  optical  phase 

{<I>  =  Q)  and  weak  field  amplitude  |a|  =  =  2  .  The  separatrix  is  plotted  in  red  and  the 

paths  of  the  electron  are  plotted  with  a  yellow  to  red  gradient.  The  separatrix  is  the  path 
that  separates  closed  orbits  (analogous  to  a  pendulum  swinging  back  and  forth)  from 
open  orbits  (analogous  to  the  pendulum  going  over  the  top  and  swinging  in  a  full  circle) 
and  is  given  by 

v"  =  2|fl|[l  +  sin(4'+^^)].  (III.24) 

The  height  of  the  separatrix  is  determined  by  the  strength  of  the  optical  field 
=  2V2  .  The  final  electron  positions  are  represented  in  blue.  This  plot  illustrates 

there  is  an  optimal  energy  exchange  at  resonance,  with  some  electrons  gaining  and  some 
losing  energy.  The  higher  energy  electrons  catch  up  to  those  with  lower  energy,  and 
bunching  occurs  as  shown  by  the  decrease  of  spacing  between  the  electrons  toward  the 
center  in  the  figure  below.  However,  in  this  case,  an  equal  number  of  electrons  gain  and 
lose  energy,  resulting  in  a  net  energy  change  in  the  electron  beam  of  nearly  zero.  There 
is  then  negligible  gain  G  in  the  intensity  of  the  laser  light. 
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***  FBL  Phase  S^ace  JSvoIution  *** 
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Figure  6.  FEL  Phase  Space  Plot  at  Resonance  (vo  =  0)  ([From  Ref.  [4].) 


A  more  desirable  phase  space  plot  of  an  operating  FEL  is  shown  in  Figure  7.  In 
this  simulation,  1000  sample  electrons  start  out  at  energies  slightly  above  resonance 
( Vp  =  3 ).  There  is  energy  exchange  between  the  electrons  and  the  light,  but  more 

electrons  lose  energy  than  gain  energy,  which  is  desirable  since  the  energy  lost  by  the 
electrons  is  gained  by  the  light.  There  is  also  a  bunching  of  electrons  that  are  losing 
energy,  leading  to  coherent  radiation,  another  desired  characteristic  of  a  free  electron 
laser. 


In  (l-i-G) 


0.15 


0.05 


Figure  7. 


Phase  Space  Above  Resonance  (vo  =  3)  ([From  Ref  [4].) 
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G.  DIFFRACTION  OF  THE  LASER  BEAM 


Propagation  of  a  light  wave  can  be  mathematically  described  using  the  wave 
equation  given  by 


f 


—2 

V 


V 


c'  de  ^ 


E{x,t)  =  0 . 


(111.25) 


This  equation  involves  four  second  order  derivatives  in  (x,  t)  for  the  propagation  of  the 
light,  with  E{x,t)  representing  the  electric  field  of  the  light.  In  a  typical  operating  FEL, 
the  light  is  highly  coherent  and  the  electric  field  due  to  the  light  is  slowly  varying  in  the 
direction  of  travel  (z  axis)  over  an  optical  wavelength  and  slowly  varying  in  time  over  the 
optical  frequency.  The  complex  electric  field  is  then  represented  by 

E{x,t)  =  a{x,t)e'"  =  \a{x,t)\e'*^^^'‘^e'"  [7].  (III.26) 


In  this  equation,  a  =  kz- cot  represents  the  carrier  wave  and  the  optical  field  is 
represented  bya  =  If  the  optical  field  is  constant  in  magnitude,  then  the  electric 

field  E  describes  a  simple  plane  wave  traveling  in  the  z  direction  containing  a  single 
frequency  co .  Because  the  laser  amplitude  and  phase  are  slowly  varying,  we  can  make 
the  following  assumptions  describing  the  two 

d\a\  I  I  d\a\  ,  \  \  dc/)  ,  dcj)  ,  , 
dt  dz  dt  dz 

Substituting  the  expression  for  E(x,t)  into  the  wave  equation,  gives 


VU  + 


d^a 


da 


A 


,  +2ik  —  -k^a 
dz 


d^a 


da 


A 


■-2ico—-co^a 
df  dt 


=0, 


(III.27) 


where  ^  +  ^  is  the  transverse  Laplacian.  Because  of  the  slowly  varying 

dx  dy 

amplitude  and  phase  assumptions,  we  may  simplify  the  equation  by  neglecting  the  second 
derivative  terms.  Multiplying  the  simplified  equation  by  the  factor  e“'“  results  in 
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V\+2ik 


^  5  1  5  ^ 

ydz  c  dt j 


a(x,t)  =  0. 


(111.28) 


This  equation  is  simplified  further  by  introducing  another  coordinate  u  =  z-ct  that 
follows  the  wavefront.  The  partial  derivatives  are  now 


d  _du  d  dt  d  _  d 
dz  dz  du  dz  dt  du  ’ 
\  d  _  \  f du  d  ^dt  d'^ 
c  dt  c\dt  du  dt  dt  y 


d  ^  \  d 
du  c  dt 


Substitution  of  the  derivative  forms  using  the  new  coordinate  now  gives  the  wave 
equation  a  new  form 


vi  + 


2ik  d 
c  dt 


a(x,t)  =  0 . 


(III.29) 


A  dimensionless  time  variable  is  then  defined  relating  to  the  travel  time  of  the 
laser  light  along  the  undulator  T  =  ct/ L,  where  L  is  the  undulator  length.  This  time 
variable  ranges  from  r  =  0  at  the  beginning  of  the  undulator  to  r  =  1  at  the  end  of  the 
undulator.  Introduction  of  this  variable  results  in 


Y-zr) 

1^2  5 

v,  + — 

.1  2k  ) 

’  dr 

a{x,t)  =  0 . 


(III.30) 


The  first  term  in  equation  III. 30  is  the  diffraction  term,  and  its  magnitude  depends 
on  the  transverse  area  of  the  laser  beam,  its  wavelength,  X  =  2n !  k  and  the  range  to  the 
target.  The  Rayleigh  length  is  defined  as  the  distance  from  the  mode  waist  to  where 
the  area  of  the  optical  mode  doubles  due  to  diffraction.  This  length  relates  the  optical 
mode  waist  radius  to  the  wavelength  of  the  light  by  the  following:  I Z.  In 

defining  dimensionless  variables  for  transverse  coordinates,  is  taken  to  be  equal  to  the 
length  of  the  undulator  L.  The  characteristic  mode  waist  radius  is  ^/Zz7^  =  yj2Llk  . 
When  the  optical  mode  radius  is  much  larger  than  \IlZI  tt  ,  then  diffraction  has  a  small 

effect  over  the  range  of  propagation,  if  it  is  smaller  than  -sJlZI  tt  ,  diffraction  has  a  big 
effect  in  propagation  to  the  target.  Dimensionless  coordinates  based  on  the  characteristic 
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mode  waist  radius  are  x  =  xylk/lL  and  y  =  yyjk/lL  .  Substitution  of  these 
dimensionless  variables  yields  a  totally  dimensionless  form  of  the  wave  equation 


da 

dr 


(111.31) 


d^  d^ 

where  - ^ .  When  the  transverse  eomponents  (x,y )  of  the  laser  beam  are  on 

dx  dy 

the  order  of  unity  or  smaller,  diffraction  has  a  significant  effect  over  the  propagation 
range  ( r  =  0  — >  1 );  if  they  are  much  greater  than  unity,  diffraction  has  little  effect  and  can 
be  neglected.  Diffraction  causes  both  the  amplitude  and  phase  of  the  optical  field  to 
evolve. 


H.  FEL  WAVE  EQUATION 

The  FEL  wave  equation  describes  the  evolution  of  the  optical  field  interacting 
with  the  electron  beam  within  the  undulator.  The  wave  equation  with  a  current  source 

J_^  is 


r 


V 


2 


V 


dt^  ^ 


A(x,  t)  = - (x,  t)  . 

c 


(III.32) 


In  this  equation,  A(x,  t)  represents  the  optical  vector  potential.  The  laser  electric  field  can 
be  obtained  from  this  potential  by 

1  dA 

E  =  — — .  (III.33) 

c  dt 

As  before,  the  optical  field  is  coherent  laser  light  so  E{x,t)  is  complex  and  slowly 
varying  in  amplitude  and  phase  over  a  laser  wavelength  and  slowly  varying  in  time 
compared  to  the  optical  frequency.  The  vector  potential  is  given  by 

A{x,  t)  =  .  (III.34) 

k 

As  before,  a  =  kz-cot  represents  the  carrier  wave,  and  s  is  the  laser  polarization  vector 
(circularly  polarized  due  to  the  helical  undulator,  s  =  (-/,  1, 0) ).  Using  the  slowly  varying 
amplitude  and  phase  assumptions  the  complete  wave  equation  becomes 
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The  current  source  ,  where  Vj^  =  c/?j^  is  due  to  the  transverse  motion 

of  the  electrons  in  the  undulator.  Introducing  dimensionless  variables  for  time  and 
transverse  coordinates  (t  =  ct  I  L ,  x  =  x(k / 2L) ,  and  y  =  y(k / 2L) )  and  multiplying  by 

{-ATrNeKJ}  /  )  gives  the  dimensionless  wave  equation 

-jVi+A  <,(*,/)  =  -< (III.37) 

In  this  equation,  a  =  \a\e''^  represents  the  dimensionless  optical  field  amplitude 

|a|  =  4;7-AAeiS|£'|  /  ,  j  is  the  dimensionless  current  j  =  p!  ,  p 

represents  the  density  of  electrons  in  a  small  volume  element,  and  <  >  represents  the 

average  over  all  sample  electrons  in  a  small  volume  element.  If  there  is  little  diffraction 
in  the  undulator  and  the  electron  and  optical  beams  overlap,  the  wave  equation  simplifies 
to 

0 

a  =  -j<  e-'^  >  .  (111.38) 

When  TT ,  the  PEL  is  operating  in  the  high-gain  regime;  when  j  <7i ,  the  FEE  is  in 

the  low-gain  regime.  The  factor  indicates  the  amount  of  electron  bunching  in 

the  undulator. 
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IV.  FREE  ELECTRON  LASER  ADVANCED  TOPICS 


A.  WEAK  FIELD  GAIN  THEORY 

Gain  in  the  free  electron  laser  is  achieved  through  the  interaction  of  the  bunched 
electron  beam  with  the  optical  field  in  the  undulator.  This  section  will  derive  gain  in 
weak  optical  fields  (|a|<iC;T).  The  following  expressions  are  valid  for  weak  or  strong 
optical  fields 


d  a 
dr 


=  -y<cos(4'  +  ^^)>, 


(IV.l) 


^  =  f4<sin(^  +  <zi)>,  (IV.2) 

rl 

where  <  >  represents  the  average  of  all  the  electrons  within  an  optical  wavelength  of  the 
beam,  and  describe  the  change  in  optical  amplitude  and  phase  during  PEL  operation  due 
to  the  electron  beam  interaction.  In  examining  these  equations,  optimal  gain  occurs  when 
the  factor  <  cos(4’  +  </>)>  is  approximately  negative  one.  This  condition  occurs  when  the 
electrons  are  bunched  around  the  phase  +  (I)kk  . 

There  are  two  distinct  gain  regimes,  low  and  high.  These  regimes  are 
characterized  by  the  dimensionless  current  density  j,  describing  the  coupling  between  the 
electron  beam  and  the  optical  field,  as  shown  in  the  PEL  wave  equation  (equation  III.42) 
and  the  above  equations.  If  j  <7i ,  the  system  is  in  the  low  gain  regime;  if  7  »  ,  it  is  in 

the  high  gain  regime. 

1.  Low  Gain  Regime 

In  the  low  gain  regime,  j  is  small  and  there  is  not  an  appreciable  change  in  the 
optical  field  amplitude  or  phase  as  the  electrons  pass  through  the  undulator.  This  results 
in  an  approximately  fixed  separatrix  when  the  electron  evolution  is  plotted  in  phase 
space.  The  small  change  of  the  optical  amplitude  and  phase  as  the  electron  beam  goes 
through  the  undulator  (r  =  0— >^l)is  given  by 

\a(T)\  = 


1  2- 

i"^l- 


■  2 cos( Vq t)-  VqT sin( t)) 

J 


(IV.3) 
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(IV.4) 


(K'^)  =  j 


r2sin(vQr) 


-  v,r(l  +  cos(Vor)^ 

o  3  L'J’ 

2^0  ) 


where,  a,,  represents  the  initial  optical  field  amplitude  and  is  the  initial  phase  velocity 
of  the  electron  beam.  These  equations  are  obtained  by  expanding  and  v  in  powers  of 
\a\  =  Uq  .  Gain  in  the  PEL  system  is  sensitive  to  initial  phase  velocity  or  energy  of  the 
electron  beam.  The  electron  beam  starting  at  resonance  ( =  0 ,  see  Figure  6)  yields 
negligible  gain,  and  initial  phase  velocities  just  below  resonance  (v,,  <  0)  yield  negative 

gain  (absorption)  in  the  system.  The  resulting  gain  as  the  electrons  go  through  the 
undulator  is  given  by 


2  -  2 cos(Vor)  -  Vorsin(Vor)j 


(IV.5) 


The  weak-field,  low  gain  and  phase  spectra  as  functions  of  initial  phase  velocity 
are  shown  below  in  Figure  8.  The  gain  spectrum  is  anti-symmetric  about  the  resonance 
frequency  with  positive  gain  for  phase  velocities  just  above  resonance  (see  Figure  7). 
The  maximum  gain  of  about  14  percent  is  achieved  with  an  initial  phase  velocity  of 
approximately  tt  with  7  =  1.  The  phase  spectrum  is  symmetric  about  resonance,  where 
the  optical  phase  shift  is  at  a  maximum. 

+  *  +  *  Caiii  aad  Phase  Curves 
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Gain  and  Phase  Spectra  for  Low-Gain  ([From  Ref  [4].) 
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Figure  8. 


2.  High  Gain  Regime 

Unlike  the  low  gain  regime,  when  the  dimensionless  current  density  is  large 
( 7  »  ),  it  strongly  affects  the  optical  field  amplitude  and  phase  as  the  electron  beam 

goes  through  the  undulator.  The  optical  field  amplitude  grows  exponentially,  and  the 
phase  grows  linearly.  The  optical  field  amplitude  and  phase  during  the  electron  beam 
interaction  with  the  optical  field  are  given  by 


=  ’  ,  and 

V  3 

(IV.6) 

^{t)={J/2)''\/2  [7], 

(IV.7) 

Because  of  the  large  changes  in  the  optical  field  phase,  the  separatrix  in  phase 
space  is  shifted  and  the  optimum  bunching  location  changes.  The  gain  spectrum  changes 
as  well  and  the  final  gain  is  not  as  sensitive  to  initial  phase  velocity  as  was  the  case  in  the 
low  gain  regime.  A  simulation  of  an  PEL  operating  in  the  high  gain  regime  (j  =  200 ) 
with  weak  fields  is  shown  in  Figure  9.  In  comparison  to  Figure  7,  where  the  bunching  is 
at  a  phase  of  ^  ,  the  separatrix  and  bunching  location  has  shifted  to  the  left  by  about 

TT .  In  this  simulation,  the  electron  beam  started  at  resonance  ( =  0 ),  and  yet  there  was 
considerable  gain. 

FEJL  Ehase  space  E-Jolution  *** 
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Phase  Space  of  High  Gain  Regime  ([From  Ref  [4].) 
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Figure  9. 


The  expression  for  gain  as  the  eleetron  beam  goes  through  the  undulator  in  the 
high-gain  regime  is  given  by 

=  [7],  (IV.8) 

9 

The  gain  and  phase  spectra  of  the  high-gain  regime  are  shown  in  Figure  10,  where 
the  gain  is  almost  symmetric  about  the  resonance  phase  velocity.  Discontinuities  in  the 
phase  spectrum  are  due  to  periodic  boundary  conditions  {-7i<<j)<K^. 


Giiii  aad  Phase  Curves 
j-20D  SI-30 
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Figure  10.  Gain  and  Phase  Spectra  for  High  Gain  ([From  Ref  [4].) 

B.  OPTICAL  GUIDING 

During  FEL  operation,  the  electron  beam  interaction  opposes  the  normal 
diffraction  effects  of  light  propagation  within  the  undulator.  The  optical  mode  is 
typically  wider  than  the  electron  beam  that  interacts  with  it.  The  electrons  amplify  the 
area  of  the  optical  mode  that  they  interact  with.  Therefore,  the  edges  of  the  optical  mode 
are  not  amplified,  and  as  the  FEL  interaction  continues,  the  central  part  of  the  mode  that 
is  interacting  with  the  beam  grows  in  amplitude.  This  interaction  results  in  a  focusing 
effect  on  the  optical  field  in  the  undulator  and  counteracts  the  normal  diffraction  of  the 
optical  field.  This  is  known  as  gain  guiding.  In  addition,  with  increasing  gain,  there  is 
increasing  optical  field  phase  as  shown  in  equation  IV.7.  This  optical  field  phase  change 
counteracts  the  phase  change  that  would  occur  due  to  diffraction.  This  is  known  as 
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optical  phase  guiding.  Both  optical  phase  and  gain  guiding  benefit  the  PEL  interaction, 
preventing  diffraction  and  focusing  the  optical  field  so  it  overlaps  the  electron  beam 
throughout  the  undulator.  This  feature  allows  the  use  of  a  long  undulator,  a  vital 
component  of  a  high-gain  amplifier  PEL. 

C.  BETATRON  MOTION 

In  addition  to  optical  guiding,  there  is  a  focusing  of  the  electron  beam  itself  within 
the  undulator.  In  the  basic  theory  chapter,  the  electrons  were  assumed  to  be  perfectly 
injected  into  the  undulator.  In  the  case  of  imperfect  injection,  the  electrons  have  variable 
displacements  and  angles  from  the  longitudinal  axis  of  the  undulator.  A  displacement 
from  the  axis  results  in  an  imbalance  of  forces  due  to  the  undulator  magnets,  resulting  in 
a  low  frequency  sinusoidal  motion  known  as  betatron  oscillation.  The  dimensionless 
betatron  oscillation  frequency  is  given  by 


IttNK 

C0p= -  [7]. 

r 


(IV.9) 


In  typical  PELs,  the  frequency  is  on  the  order  of  Ik  so  there  is  about  one  betatron 
oscillation  over  many  periods  of  the  undulator  as  r  =  0  ^  1 .  An  illustration  of  betatron 
oscillation  is  shown  in  Pigure  11. 


Pigure  1 1 .  Betatron  Pocusing  and  Oscillation 


In  this  figure,  the  electron  beam  (shown  in  red)  enters  at  an  angle  that  would  lead 

to  transverse  beam  drift  along  the  undulator.  However,  the  magnetic  fields  cause  the 

beam  to  arch  downwards  back  towards  the  center  of  the  undulator.  The  beam’s  final 

position  is  slightly  below  the  center  of  the  undulator.  If  the  undulator  were  longer,  one 

would  observe  the  beam  bending  upward  back  toward  the  center.  As  shown  in  the  figure, 

the  frequency  of  oscillation  is  low  and  there  is  typically  less  than  one  oscillation  within 
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the  undulator  length.  The  amplitude  of  this  oscillation  for  each  electron  depends  on  the 
beam’s  initial  displacement  and  angle  from  the  axis  of  the  undulator.  This  betatron 
oscillation  is  vital  in  keeping  the  beam  focused  within  the  undulator,  and  allowing  for 
imperfect  electron  beam  injection. 

D.  STRONG  FIELD  GAIN  THEORY 

The  initial  strength  of  the  optical  field  also  plays  an  important  role  in  the  final 
gain  of  the  system.  In  both  the  low  and  high  gain  regimes,  a  strong  optical  field  makes 
the  electron  beam  overbunch  during  travel  through  the  undulator.  This  process  is  called 
saturation.  Overbunching,  as  the  electrons  travel  further  through  the  undulator,  results  in 
energy  being  taken  away  from  the  optical  field,  reducing  overall  gain  of  the  system.  In 
the  low-gain  regime,  saturation  is  achieved  at  much  lower  initial  field  amplitudes 
(a^  than  the  high-gain  regime  (a^  k2(J/2Y^^)  [7].  An  example  of  overbunching 

(saturation)  in  the  low  gain  regime  is  shown  in  Figure  II.  When  there  is  no  saturation, 
the  phase-space  evolution  of  the  electrons  is  slower  and  the  bunch  of  electrons  ends  up  in 
a  downward  trend  bunched  at  a  phase  of  ^  as  shown  in  Figure  7.  However,  in  this 
situation,  the  optical  field  is  stronger,  causing  a  greater  force  on  the  electrons  and  thus 
more  evolution  in  phase  space.  The  bunch  forms  and  continues  along  closed  orbits  until 
it  starts  an  upward  trend,  taking  energy  back  from  the  optical  field.  This  saturation  can 
be  extended,  however,  to  achieve  greater  gain  by  tapering  the  undulator. 
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Low  Gain  Strong  Field  Saturation  ([From  Ref.  [4].) 
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Figure  12. 


E.  TAPERING 


To  increase  extraction  levels  in  the  amplifier  configuration  of  the  PEL,  some 
designs  have  utilized  the  concept  of  tapering  the  undulator.  Tapering  alters  the  resonance 
condition  (equation  111.22)  as  the  electron  beam  travels  through  the  undulator,  allowing 
prolonged  electron  beam-optical  field  energy  exchange.  This  prolonged  exchange  leads 
to  higher  extraction  77 ,  or  percentage  of  energy  taken  from  the  electron  beam. 

Examining  the  resonance  condition  reveals  that  there  are  several  ways  that 
resonance  can  be  affected:  changing  the  undulator  parameter,  undulator  period,  or 
electron  beam  energy.  Almost  all  of  the  tapered  designs  of  FELs  utilize  the  easiest  of  the 
three  methods:  varying  the  undulator  parameter  by  changing  the  magnetic  fields  along  the 
undulator.  This  is  typically  done  by  changing  the  separation  distance  (gap)  between  the 
magnets,  either  linearly  or  in  a  “step”  change  [8]. 

In  the  linear  taper,  a  reduced  (negative  taper)  or  increased  (positive  taper)  gap 
between  facing  magnets  is  introduced  after  the  first  period  of  the  undulator.  With  each 
undulator  period,  the  gap  is  either  reduced  or  increased  linearly  through  the  rest  of  the 
undulator.  The  increase  (or  decrease)  in  the  gap  is  very  small  (on  the  order  of  tenths  of 
millimeters).  This  technique  effectively  exposes  the  electron  beam  to  a  decreasing 
magnetic  field  as  shown  in  Figure  12  in  the  case  of  positive  linear  tapering,  and  an 
increasing  magnetic  field  as  shown  in  Figure  13  for  negative  linear  tapering. 


Figure  13.  Positive  Linear  Tapered  Undulator  Field  ([From  Ref  [4].) 
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Figure  14. 
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Negative  Linear  Tapered  Undulator  Field  ([From  Ref.  [4].) 


The  step  taper  is  implemented  by  changing  the  undulator  gap  at  a  certain  point 
within  the  undulator  and  maintaining  the  reduced  or  increased  gap  through  the  rest  of  the 
undulator.  Examples  of  the  positive  and  negative  step-tapered  undulator  are  shown  in 
Figure  14  and  Figure  15. 
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Figure  16. 
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Positive  Step-Tapered  Undulator  Field  ([From  Ref.  [9].) 
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Negative  Step-Tapered  Undulator  Field  ([From  Ref  [9].) 


Tapering  changes  the  pendulum  equation  to  the  following  form, 


^  =  5  +  \a\co^{^  +  (t>)  [7], 


(IV.  10) 


S  =  -47rN- 


AK 
l  +  K^  K 
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where  S  is  the  phase  acceleration  introduced  by  the  tapering  of  the  undulator  and 
AK  /  K  represents  the  fractional  change  in  the  undulator  parameter.  Recall  that  K  is 
proportional  to  B,  the  magnetic  field  strength. 

This  phase  acceleration  plays  a  big  part  in  the  additional  extraction  that  a  system 
can  gain  due  to  the  tapering  of  the  undulator.  The  optimum  extraction  from  a  normally 
configured  undulator  is  approximately 


7 


(IV.ll) 


while  the  amount  of  extraction  from  a  tapered  undulator  is  given  by 


8;rN 


m. 


(IV.  12). 


Equation  IV.ll  shows  that  the  more  taper  increases,  the  greater  the  extraction 
levels  will  be;  however,  there  is  an  upper  limit  to  the  amount  of  tapering  that  can  be  used. 
The  upper  limit  is  set  by  the  fact  that  the  pendulum  equation  has  no  solution  when  the 
taper  phase  acceleration  is  greater  than  the  optical  field  amplitude.  Also,  a  lower  limit  on 
tapering  is  determined  by  the  maximum  deceleration  a  trapped  electron  can  undergo  in 

I  |l/2 

the  non-tapered  case.  This  is  given  by  doubling  the  separatrix  height,  2 \a\  .  The  two 

conditions  set  the  optimal  range  for  tapering  given  by 

4yfa<S<\a\  (IV.  13) 


In  phase  space,  the  separatrix  that  results  from  this  altered  pendulum  equation  is 
oriented  differently  depending  on  the  type  (positive  or  negative)  of  tapering.  Figure  17 
illustrates  the  phase  space  evolution  for  a  positively  tapered  undulator,  while  Figure  18 
depicts  the  results  of  negative  tapering. 
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*44  PEX  Phase  Space  Evolution  *** 
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Figure  17.  Positive  Linear  Taper  Phase  Spaee  Plot  (From  Ref.  [4].) 


***  TEL  Phase  Space  Evolution  *** 
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Figure  18.  Negative  Linear  taper  Phase  Space  Plot  ([From  Ref.  [4].) 

In  the  case  of  positive  taper,  some  of  the  electrons  are  trapped  inside  the 
separatrix  while  the  others  are  in  open  orbits,  while  in  the  negative  tapering  case,  none  of 
the  electrons  are  trapped  inside  the  separatrix,  but  must  go  around  it.  This  motion  around 
the  separatrix  causes  bunching  and  good  extraction.  However,  negative  taper  only  works 
well  when  the  electron  beam  starts  above  resonance. 

The  separatrix  in  phase  space  decreases  in  size  as  the  phase  acceleration  (or 


deceleration)  approaches  the  value  of  the  optical  field  amplitude.  When  the  phase 
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acceleration  is  greater,  the  separatrix  vanishes.  Without  the  separatrix,  there  are  no 
trapped  electrons  and  no  significant  bunching  occurs. 

Tapering,  when  applied  optimally  within  the  valid  range,  can  yield  significant 
improvement  in  extraction  levels  over  a  conventional  PEL,  especially  in  the  high-gain 
and  strong  field  regimes. 
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V.  BROOKHAVEN  NATIONAL  LAB  100-KW  FEL 


A.  INTRODUCTION 

The  free  electron  laser  shows  great  promise  as  a  naval  weapon,  but  is  not  yet 
mature  enough  in  its  design  to  produce  the  power  output  levels  required  for  use  in 
weapons  application.  The  most  powerful  FEL  in  the  world  today  (an  oscillator 
configuration)  is  located  at  Jefferson  Lab.  It  produces  a  continuous  wave  output  of  10 
kW,  two  orders  of  magnitude  short  of  the  MW  level  needed.  Two  competing  designs  for 
a  100  kW  level  FEL,  proposed  by  Brookhaven  and  Los  Alamos  National  Labs,  both 
amplifiers,  provide  another  vital  step  in  the  evolution  of  the  weapons  level  design.  This 
chapter  will  focus  on  the  Brookhaven  design,  which  utilizes  a  conventional,  non-tapered 
undulator.  The  parameters  of  the  design  are  shown  in  Table  1. 


Frequency  (MHz) 

704 

Wiggler  Input  Energy  Spread  (%) 

0.02 

Bunch  Charge  (nC) 

1.4 

Wiggler  Period  (mm) 

32.5 

PRF  (MHz) 

357 

Wiggler  Length  (m) 

3.9 

Average  Beam  Current  (A) 

0.5 

Number  of  Wiggler  Periods 

120 

Injector  Energy  (MeV) 

2.5 

Wiggler  Gap  (mm) 

20 

Accelerating  Gradient  (MV/m) 

20 

Peak  Beam  Current  (A) 

500 

Transverse  rms  Emittance  (qm) 

5.0 

Extraction  (%) 

0.5 

Longitudinal  rms  Emittance  (keV-psec) 

38 

IR  Output  Power  (kW  @  1  pm) 

200 

Wiggler  Energy  (MeV) 

80 

Bunch  Length  (psec) 

2.8 

Wiggler  Beam  Power  (MW) 

40 

Beam  Dump  Power  (MW) 

1.25 

le  1 .  Brookhaven  National  Lab  (E 

iNL)] 

Design  Parameters  ([After  Ref  | 

[10].) 

B.  PERFORMANCE  OF  THE  BROOKHAVEN  FEL  DESIGN 

Computer  simulations  using  codes  developed  at  the  Naval  Postgraduate  School 
were  run  to  determine  the  performance  of  the  design  with  the  electron  beam  misaligned 
in  various  ways. 
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1.  Weak-Field  Gain  Spectrum 

The  Brookhaven  PEL  design  was  first  simulated  with  an  aligned  electron  beam 
injection  to  study  its  gain  and  extraction  spectrum.  To  obtain  the  gain  spectrum,  the 
electron  beam  phase  velocity  vo  (proportional  to  energy)  was  varied  and  the  optical  field 
kept  in  the  weak  field  regime  ( {j  /  2)^^^ ).  The  gain  spectrum  is  shown  in  Figure  19. 


Gain  vs.  v 
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Figure  19.  BNL  FEE  Weak-Field  Gain  Spectrum 


Examining  the  gain  spectrum,  one  sees  the  broad  spectrum  indicative  of  a  high  gain  FEE, 
with  a  maximum  of  G  ~ 12000  at  Vq=4.  This  phase  velocity  corresponds  to  a  optical 

wavelength  shift  of  0.5%  above  resonance.  The  effects  of  high  gain  optical  guiding  can 
be  examined  at  the  peak  of  the  gain  spectrum.  A  picture  of  the  interacting  electron  beam 
(shown  in  red)  and  optical  field  (shown  in  yellow)  within  the  undulator  is  shown  in 
Figure  20. 
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Figure  20.  Simulation  of  Gain  Guiding  at  Peak  Gain  of  BNL  Design 

At  the  beginning  of  the  undulator,  where  gain  is  low,  the  initial  optical  mode  is 
fairly  wide.  As  the  optical  mode  proceeds  through  the  undulator  and  interacts  with  the 
electron  beam,  gain  guiding  amplifies  the  central  part  of  the  mode  and  counteracts  the 
phase  change  due  to  diffraction,  focusing  the  mode  in  the  center  of  the  undulator.  The 
resulting  amplified  mode  is  significantly  narrower  at  the  end  of  the  undulator.  This 
reduced  the  optical  mode  size,  so  that  it  will  diffract  to  a  larger  area  at  the  first 
subsequent  optic,  resulting  in  a  smaller  intensity  on  the  optic.  The  minimization  of 
intensity  on  optics  is  an  important  issue  in  the  design  of  any  high-power  FEL. 

2.  Extraction  Spectrum 

Simulations  at  various  initial  electron  beam  phase  velocities  vo  determined  the 
maximum  extraction  rj ,  defined  as  the  output  laser  power  divided  by  the  power  in  the 
electron  beam  entering  the  undulator.  The  simulations  result  in  the  extraction  spectrum 
of  Figure  21. 
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Extraction  vs.  v 
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The  maximum  extraction  that  this  design  will  achieve  is  approximately  one 
percent.  This  maximum  occurs  at  an  electron  beam  initial  phase  velocity  of  =  1 1 , 

corresponding  to  a  wavelength  shift  of  1.5%  away  from  resonance.  With  an  electron 
beam  power  of  40  MW,  this  extraction  yields  an  output  power  of  400  kW,  which  is  four 
times  the  design  specification  of  the  PEL.  A  minimum  extraction  of  only  0.25%  is 
required  to  yield  the  design  specified  100  kW  output  power. 

From  the  optimal  extraction  simulation,  the  electron  beam-optical  field  interaction 
yields  the  phase  space  picture  of  Figure  22,  where  there  is  an  induced  electron  beam 
energy  spread  of  4%.  The  vertical  axis  is  proportional  to  the  electron  energy  while  the 
horizontal  axis  is  the  electron  phase  within  an  optical  wavelength. 


-n/2  Z  37t/2 

Figure  22.  BNL  FEE  Optimal  Phase  Space 
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Gain  is  not  optimal  when  extraction  is  maximized,  but  the  PEL  still  reaps  the 
benefits  of  gain  guiding  as  shown  in  Figure  23.  The  Gaussian  shape  in  red  at  the  right 
end  of  the  undulator  represents  the  optical  mode  size  due  to  normal  free  space  diffraction. 
The  yellow  is  the  actual  mode.  There  is  seen  to  be  significant  focusing  due  to  the 
electron  beam  interaction. 


Figure  23.  Gain  Guiding  in  Optimal  Extraction  Simulation  of  BNL  FEE. 

Optical  guiding  in  the  undulator  leads  to  increased  diffraction  and  a  bigger  mode 
radius  at  the  first  optic  as  shown  in  Figure  24.  The  theoretical  optical  mode  size  due  to 
diffraction  is  shown  in  red,  while  the  actual  mode  (approximately  twice  the  size)  is 
shown  in  yellow. 


-35  X  35 

Figure  24.  Optical  Mode  at  First  Optic  (27  m  away) 

The  spot  area  on  the  first  optic  spaced  at  27  meters  from  the  beginning  of  the 
undulator  is  approximately  12  cm^  so  the  resulting  intensity  at  the  first  optic  is  30 


39 


kW/cm^.  The  x  and  y  axes  in  Figure  24  refer  to  the  transverse  extent  of  the  optical  mode 
intensity  shown  by  the  brighter  area  in  the  background  of  the  figure. 

3.  Electron  Beam  Shift 

Vibration  tolerances  of  the  Brookhaven  FEL  design  are  analyzed  in  the  next  two 
sections.  Onboard  ship,  vibrations  affecting  the  FEL  have  an  extremely  low  frequency 
when  compared  to  the  pulse  repetition  frequency  of  the  optical  beam  in  the  undulator. 
For  this  reason,  vibration  effects  (shift  and/or  tilt  of  the  electron  beam)  can  be  modeled  as 
static  conditions  in  the  Brookhaven  FEL. 

Vibration,  in  essence,  results  in  imperfect  electron  beam  injection  into  the 
undulator.  The  case  of  a  vertical  transverse  shift  of  the  electron  beam  away  from  the 
undulator  axis  is  studied  first.  Many  simulations  optimize  the  extraction  in  vo  at  various 
amounts  of  electron  beam  shift  off  the  undulator  axis,  producing  the  results  shown  in 
Figure  25. 


Extraction  vs.  Electron  Beam  Shift 


As  expected,  performance  of  the  system  is  degraded  and  extraction  decreases  as 
the  electron  beam  shift  increases.  A  maximum  tolerable  shift  of  approximately  1.2  mm 
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would  still  allow  the  PEL  to  produce  the  specified  output  power  of  100  kW.  Figure  26 
shows  the  effect  of  a  large  electron  beam  shift  in  the  undulator. 


Figure  26.  Electron  Beam  Shift  in  the  Undulator 


Betatron  focusing  deflects  the  electron  beam  back  to  the  center  of  the  undulator. 
There  is  also  optical  guiding  as  the  beam  bends  in  the  undulator.  These  effects, 
unfortunately,  result  in  a  distorted,  non-symmetric  optical  mode  exiting  the  undulator. 
Experimental  FEE  systems  with  active  alignment  can  hold  the  electron  beam 
displacement  yo  to  within  0.05  mm,  so  there  is  expected  to  be  only  a  small  amount  of 
degradation  due  to  vibration. 

4.  Electron  Beam  Tilt 

The  Brookhaven  FEE  design  is  also  subjected  to  another  effect  of  vibration:  the 
electron  beam  entering  the  undulator  an  angle.  Simulations  were  run  (and  optimized  with 
respect  to  )  for  varying  angles  to  determine  the  effect  of  the  angle  of  the  electron  beam 

at  the  middle  of  the  undulator  on  extraction.  The  resulting  graph  is  shown  in  Figure  27. 
Each  value  of  extraction  is  evaluated  at  the  optimum  vq. 
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Extraction  vs.  Electron  Beam  Tilt  fr^  =  0.5) 

p 


The  performance  of  the  system  is  hampered  by  the  angle  of  the  electron  beam,  as 
shown  by  the  decreasing  extraction  with  increasing  angle  of  tilt.  Performance  is 
degraded  most  severely  (i.e.  steeper  slope)  with  shifts  of  approximately  0.6  mrad  or 
greater.  A  maximum  electron  beam  tilt  of  approximately  0.9  mrad  still  allows  the  system 
to  produce  designated  output  power.  A  simulation  at  the  maximum  allowable  tilt  (Figure 
28)  illustrates  this  effect  with  electron  beam  focusing  at  the  midway  point  of  the 
undulator. 


Figure  28.  Electron  Beam  Tilt  Midway  Through  Undulator 


Betatron  motion  keeps  the  electron  beam  focused  enough  in  the  undulator,  but  the 
output  optical  mode  is  shifted  vertically  upon  exiting  the  undulator  due  to  optical  guiding 
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by  the  angled  electron  beam.  These  effects  result  in  a  poorly  shaped  optical  mode  at  the 
first  optic.  An  experimental  PEL  system  with  active  alignment  is  expected  to  have  an 
angular  tolerance  of  1 0  prad,  so  that  only  a  very  small  amount  of  degradation  is  actually 
anticipated. 

The  same  procedure  is  then  performed  with  the  electron  beam  angled  upon  entry 
into  the  undulator.  The  result  of  these  simulations  is  shown  in  Figure  23,  with  extraction 
at  each  point  optimized  in  vq. 


Extraction  vs.  Electron  Beam  Tilt  =  0) 


As  seen  in  the  prior  case,  increasing  the  angle  degrades  performance,  but  the 
decrease  in  extraction  is  nearly  linear.  In  addition,  the  tolerance  of  electron  beam  tilt  is 
increased  to  approximately  1  mrad.  Figure  30  illustrates  the  effect  of  tilt  with  electron 
beam  focus  at  the  beginning  of  the  undulator. 
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Figure  30.  Electron  Beam  Tilt  at  Beginning  of  Undulator 


With  focusing  at  the  beginning  of  the  undulator,  betatron  focusing  works  better 
and  the  mode  is  not  shifted,  and  is  nearly  symmetric.  This  results  in  a  useable  optical 
mode  at  the  first  optic.  But,  again,  with  active  alignment,  such  large  angular  deflection 
would  not  occur. 

In  the  case  of  all  vibration  effects  studied,  experiments  running  today  are  able, 
through  active  alignment,  to  keep  well  within  the  thresholds,  making  vibration  a  non¬ 
issue  to  the  Brookhaven  FEE. 
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VI.  LOS  ALAMOS  NATIONAL  LAB  100  KW  FEL 


A.  INTRODUCTION 

Another  design  for  a  100  kW  FEL  has  been  proposed  by  Los  Alamos  National 
Lab  (LANL).  The  parameters  for  the  design  are  listed  in  Table  3  below.  Non-tapered 
and  linearly  tapered  variants  of  this  design  are  analyzed  in  this  chapter. 


Beam  Energy 

80.8  MeV 

Matched  Beam  Radius 

0.27  mm 

Peak  Current 

1  kA 

Beam  Radius  at  Wiggler  Entrance 

0.11  mm 

Rms  Emittance 

10  mm-mrad 

LANL  Extraction 

1.9% 

rms  Energy  Spread 

0.25 

Wiggler  Length 

2.4  m 

Wiggler  Period 

2.18  cm 

Wiggler  Parameter 

1.2 

Table  2.  Los  Alamos  National  Lab  1 

Design  Parameters  ([After  Ref  [8].) 

B.  NON-TAPERED  DESIGN  PERFORMANCE 
1.  Extraction  Spectrum 

As  a  basis  for  comparison,  a  non-tapered  version  of  the  Los  Alamos  design  is 
simulated  and  its  performance  optimized.  These  simulations  result  in  the  extraction 
spectrum  shown  in  Figure  3 1 . 


0.8 
0.7 
0.6 
0.5 

£ 

'0.4 
0.3 
0.2 
0.1 

0  5  10  15  20 

V 

o 

LANL  Non-Tapered  Extraction  Spectrum 
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Figure  31. 


The  extraction  spectrum  peaks  at  (vo=  15),  corresponding  to  a  wavelength  shift  of 
2.1%  away  from  resonance.  The  maximum  extraction  obtained  by  this  design  is 
approximately  0.75%,  producing  an  output  power  of  60  kW  from  the  8  MW  electron 
beam.  This  output  fails  to  meet  design  specifications,  showing  that  tapering  is  needed  to 
boost  performance  of  the  design. 

C.  LINEAR  TAPERED  DESIGN  PERFORMANCE 
1.  Weak-Field  Gain  Spectrum 

An  alternate  design  was  explored  using  a  linear  taper.  Several  simulations  were 
run  to  optimize  the  amount  of  taper  and  the  position  in  the  undulator  where  the  taper 
begins.  For  the  optimized  design,  the  taper  begins  at  a  point  60%  of  the  way  down  the 
undulator  length  with  a  maximum  change  AB I B  of  1 8%  in  the  undulator  magnetic  field. 
This  design  was  then  analyzed  for  performance  and  to  determine  vibration  tolerances. 
The  weak  field  gain  spectrum  is  shown  in  Figure  32. 


Gain  vs.  v 
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Figure  32.  LANL  Linear  Tapered  FEL  Weak-Field  Gain  Spectrum 

The  gain  spectrum  is  very  broad  in  vo  with  maximum  gain  of  approximately  900 
occurring  at  (vo  =  4.5),  corresponding  to  a  wavelength  shift  of  0.65%  away  from 
resonance  (vo  =  0). 
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2.  Extraction  Spectrum 

The  design  is  simulated  and  optimized  to  determine  the  maximum  extraction. 
The  extraction  spectrum  is  shown  in  Figure  33. 


Extraction  vs.  v 
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Figure  33.  LANL  Linear  Tapered  FEL  Extraction  Spectrum 

With  a  minimum  extraction  of  1.2%  required  for  an  output  of  100  kW,  this 
design’s  extraction  spectrum  peaks  at  (Vq  =  6.5),  corresponding  to  a  wavelength  shift  of 
0.94%  away  from  resonance.  The  maximum  extraction  is  found  to  be  1.66%,  resulting  in 
an  output  power  of  133  kW.  This  power  level  is  more  than  double  that  of  the  non-tapered 
design,  showing  the  benefits  of  a  tapered  undulator. 

Figure  34  shows  that  this  design,  when  operating  at  optimum,  benefits  from 
optical  guiding,  resulting  in  a  smaller  optical  mode  waist  centered  at  the  end  of  the 
undulator. 
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Figure  34.  Optical  Guidance  in  a  Linear  Tapered  LANL  FEL  Design 

As  a  benefit  of  optical  guiding,  the  output  laser  beam  diffracts  to  a  larger  size  in 
free  space  after  the  interaction  in  the  undulator,  resulting  in  a  larger  mode  waist  at  the 
first  optic  and  lower  intensity.  The  mode  waist  is  approximately  2.2  cm  at  the  first  optic 
24  meters  away,  and  the  resulting  intensity  is  approximately  9  kW/cm^.  The  resulting 
laser  mode  at  the  first  optic  is  pictured  in  Figure  35. 
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Figure  35.  Optical  Mode  at  the  First  Optical  Element  (24  m) 

This  low  intensity  level  poses  no  problems  for  the  mirrors  used  in  FEL  systems 
currently  in  operation  which  are  able  to  handle  intensities  of  100  kW/cm  for  extended 
periods  of  time. 

The  phase  space  plot  shown  in  Figure  36  differs  from  the  Brookhaven  design 
(Figure  22)  with  approximately  half  of  the  electrons  trapped  and  half  free.  The  tapered 
undulator  creates  an  increased  induced  energy  spread  of  approximately  6%. 
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Figure  36.  Phase  Space  Plot  of  Linear  Tapered  LANL  Design 

3.  Electron  Beam  Shift 

The  linear-taper  design  was  simulated  under  varying  conditions  of  electron  beam 
shift  and  tilt.  The  result  of  several  simulations  with  varying  amounts  of  electron  beam 
shift  yo  is  shown  below  in  Figure  37. 


Extraction  vs.  Electron  Beam  Shift 


Figure  37.  LANL  Linear  Tapered  FEL  Extraction  vs.  Electron  Beam  Shift 

The  EANL  design  must  extract  1 .2%  of  the  electron  beam  energy  to  achieve  the 
rated  power  output  of  100  kW.  Examining  the  figure,  this  is  accomplished  with  an 
electron  beam  shift  of  up  to  approximately  0.4  mm,  a  shift  equal  to  four  times  the  radius 
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of  the  electron  beam.  A  simulation  with  the  maximum  tolerable  electron  beam  shift 
shown  in  Figure  38  illustrates  the  effects  in  the  undulator. 


Figure  38.  Maximum  Electron  Beam  Shift  on  Linear  Tapered  LANL  FEL 

Betatron  focusing  and  optical  guiding  allows  the  FEL  to  produce  rated  output,  but 
the  laser  output  mode  is  slightly  distorted  upon  exiting  the  undulator. 

4.  Electron  Beam  Tilt 

The  linearly  tapered  design  was  simulated  with  varying  degrees  of  electron  beam 
tilt  and  focus.  The  result  of  the  simulations  is  plotted  in  Figure  39. 


Extraction  vs.  Electron  Beam  Tilt  =  0.5) 
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Figure  39.  LANL  Extraction  vs.  Electron  Beam  Tilt  at  Middle  of  Undulator 
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A  maximum  tilt  of  approximately  0.5  mrad  at  the  middle  of  the  undulator  is 
allowed  before  the  PEL  fails  to  produce  the  designed  power  output.  Figure  40  illustrates 
the  electron  beam  interaction  with  the  optical  field  in  the  undulator. 


Figure  40.  Maximum  Electron  Beam  Tilt  at  Middle  of  Undulator 

Betatron  focusing  and  optical  guiding  results  in  a  smaller  symmetric,  but  slightly  shifted 
output  laser  optical  mode  exiting  the  undulator. 

Simulations  were  repeated  with  the  electron  beam  tilted  at  the  beginning  of  the 
undulator,  resulting  in  Figure  41. 


Extraction  vs.  Electron  Beam  Tilt 


Figure  41.  LANE  FEE  Extraction  vs.  Beam  Tilt  at  Start  of  Undulator 
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With  the  electron  beam  tilted  at  the  beginning  of  the  undulator,  there  is  a  decrease 
in  tolerance  of  electron  beam  angle,  with  a  maximum  tolerable  tilt  of  0.4  mrad  for  design 
output  power.  A  simulation  snapshot  of  the  extreme  limit  shown  in  Figure  42  illustrates 
the  effect  of  the  electron  beam  tilt  on  the  optical  output  of  the  undulator. 


Figure  42.  Electron  Beam  Tilt  at  the  Beginning  of  the  Undulator 


In  this  case,  the  effects  of  the  tilt  are  not  as  apparent  in  the  undulator.  Betatron 
focusing  keeps  the  electron  beam  fairly  straight  within  the  undulator.  The  optical  guiding 
results  in  the  optical  mode  centered  in  the  undulator  and  only  slightly  distorted. 

As  with  the  case  with  the  Brookhaven  design,  the  Los  Alamos  design  can  tolerate 
vibrations  such  that  the  electron  beam,  through  active  alignment,  will  not  be  shifted  or 
tilted  enough  to  hamper  performance. 

5.  Variation  of  Electron  Beam  Focus 

In  examining  the  two  cases  of  electron  beam  tilt  at  different  locations  along  the 
undulator,  there  were  higher  extraction  levels  where  the  electron  beam  was  tilted  at  the 
beginning  of  the  undulator.  This  section  examines  the  correlation  between  system 
performance  and  electron  beam  focus  points  within  the  undulator.  Simulations  are  run 
with  perfect  electron  beam  injection,  varying  the  electron  beam  focus  point  within  the 
undulator.  The  summary  of  results  is  plotted  in  Figure  30. 
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Extraction  vs.  t. 


Figure  43.  LANL  Extraction  versus  Electron  Beam  Focus  Point 

The  extraction  of  the  FEE  is  roughly  sinusoidal  with  the  position  of  electron  beam 
focus  within  the  undulator,  with  peak  extraction  of  1.9%  at  Xp  =  0.15  and  Xp  =  0.75  down 
the  length  of  the  undulator.  At  this  new  focus  the  output  is  153  kW,  a  20  kW 
improvement  over  the  baseline  design. 


53 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


54 


VII.  CONCLUSION 


We  have  examined  both  the  Brookhaven  and  Los  Alamos  designs,  and  have 
determined  that  the  Brookhaven  design  can  produce  a  higher  output,  namely  400  kW. 
The  BNL  design  did  not  allow  for  a  tapered  undulator  due  to  limits  in  acceptable  electron 
beam  energy  spread  required  for  recirculation.  With  an  increase  in  electron  beam  energy 
and  average  current,  this  design  could  achieve  an  output  of  the  desired  megawatt  level. 

We  have  also  compared  the  non-tapered  and  linear- tapered  versions  of  the  Los 
Alamos  design,  and  find  that  the  tapered  version  has  more  than  twice  the  extraction, 
demonstrating  the  benefits  of  a  tapered  undulator.  The  linearly  tapered  design  was  able 
to  produce  an  output  of  153  kW,  while  the  non-tapered  managed  only  60  kW.  As  with 
the  Brookhaven  design,  with  increased  electron  beam  energy  and  average  current,  the 
tapered  design  is  capable  of  producing  a  megawatt  level  laser  beam. 

Vibration  effects  of  electron  beam  shift  and  tilt  were,  for  the  first  time,  studied  for 
the  amplifier  configuration  of  the  PEL.  Results  for  both  the  Brookhaven  and  Los  Alamos 
designs  prove  promising,  with  experiments  currently  in  operation  able  to  hold  electron 
beam  tolerances  for  both  shift  and  tilt  to  within  an  order  of  magnitude  less  than  the 
physical  tolerances  of  the  designs  in  achieving  rated  output.  The  designs,  when 
simulated  with  shift  and  tilt  on  the  order  of  the  experimental  tolerances,  had  negligible 
performance  degradation. 

Future  work  in  this  area  includes  further  analysis  of  the  Brookhaven  FEL  to 
minimize  the  energy  spread  induced  in  the  electron  beam  during  operation  while 
maintaining  the  designed  output  level.  There  is  an  alternate  100  kW  FEL  operates  at 
Brookhaven  National  Lab  which  could  be  studied  by  us  once  we  receive  the  physical 
parameters  of  the  design. 

The  actual  design  proposed  by  Los  Alamos  National  Lab  was  a  step-tapered 
undulator.  This  design  could  be  studied  more  closely  for  optimization  and  tolerance 
effects  similar  to  those  performed  for  this  thesis. 
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Once  a  design  is  agreed  upon  that  will  yield  the  power  necessary  for  a  possible 
weapon  system,  there  are  several  issues  that  must  be  addressed  for  using  the  PEL 
onboard  a  ship.  These  include  beam  transport,  optical  element  damage  (namely  the  first 
optical  element),  and  propagation  of  the  high  energy  laser  through  the  atmosphere  to  the 
target. 
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